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ABSTRACT 

We analyze 3 epochs of ultraviolet (UV), optical and near-infrared (NIR) observa¬ 
tions of the Taurus transitional disk GM Aur using the Hubble Space Telescope Imaging 
Spectrograph (STIS) and the Infrared Telescope Facility SpeX spectrograph. Obser¬ 
vations were separated by one week and 3 months in order to study variability over 
multiple timescales. We calculate accretion rates for each epoch of observations using 
the STIS spectra and find that those separated by one week had similar accretion rates 
(~ 1 x ICC 8 Mgt/r -1 ) while the epoch obtained 3 months later had a substantially 
lower accretion rate (~ 4 x 10 -9 Mq yr~ l ). We find that the decline in accretion rate is 
caused by lower densities of material in the accretion flows, as opposed to a lower surface 
coverage of the accretion columns. During the low accretion rate epoch we also observe 
lower fluxes at both far UV (FUV) and IR wavelengths, which trace molecular gas and 
dust in the disk, respectively. We find that this can be explained by a lower dust and 
gas mass in the inner disk. We attribute the observed variability to inhomogeneities in 
the inner disk, near the corotation radius, where gas and dust may co-exist near the 
footprints of the magnetospheric flows. These FUV-NIR data offer a new perspective 
on the structure of the inner disk, the stellar magnetosphere, and their interaction. 

Subject headings: Accretion, accretion disks, Stars: Circumstellar Matter, Stars: Pre 
Main Sequence 


1. Introduction 

The early picture of static, axisymmetric circumstellar disks irradiated by a constant luminos¬ 
ity T Tauri star has evolved to include disk warps, fluctuating accretion of material in the disk 
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and ever changing emission properties of the central star, among other dynamic processes. Each 
of these manifest themsel ves in T Tauri light curves as variability . Warps block light from the 


star, causing optical dips 

(Bouvier et al. 

1999: Alencar et al. 

2010; Morales-Calderon et al. 2011; 

Codv et al. 20141 while variable accretion may be seen 

as bursts in optical and ultraviolet light 
ng widths and intensities of emission lines 

curves (Herbst et al. 1994: 

Stauffer et al. 

2014) or changi 

produced by accretion (Costigan et al. 2012; Chou et al. 

2013; 

Dupree et al. 2014). Activity cycles 

in the star may cause long timescale X-ra 
X-ray bursts are produced by solar-like fh 

y variability (A 

acela & Marino 20031), while short, strong 

ires (Wolk et a] 

. 2005). Depending on the cause, variabil- 


ity is seen as periodic or stochastic, lasting for minutes or slowly varying over a period of years. 
With all these characteristics of T Tauri light curves, often with multiple phenomena occurring 
simultaneously, explaining variability is complicated, especially if the time sampling is low or the 
wavelength coverage is limited. 


In addition, the radial structure of dust in the disk is not always continuous. Along with 


known to exist (Strom et al. 1989: 

Skrutskie et al. 1990; 

— r -* - - j. ,/ ^ 

Galvet et al. 2002: Espaillat et al. 200 o; 

Andrews et al. 20'. 

1) all thought to eventually end up as debris disks with little gas and processed 

dust ( 

Zuckerman 

2001: Chen et al. 

2006; Matthews et al. 

2014). Planet formation, photoevapo- 


ration of the disk by the central star or visco us evolution, and sett l ing have all been invoked to 


expla in these structures observed in the dust (jEspaillat et al.112014 : IZhu et al. 


2011 


Owen et al 


2012). Gaps and holes in the disk can alter interpretations of T Tauri variability, as processes rele¬ 
vant to full disks appear differently in transitional disks. For example, a warp in the disk wall at the 
sublimation temperature may not be a contributor to the light curves of transitional disks, which 
have walls at much larger radii. Therefore, varied dust distributions add an additional complication 
to the interpretation of T Tauri variability. 

Here, we analyze variability in the transitional disk GM Aur over 3 epochs of observation with 
full wavelength coverage from the far-ultraviolet (FUV) to near-infrared (NIR), all obtained within 
one day. While photometric variability studies focusing on one or multiple wavelength regimes are 
significant in number, there are few, if any, variability studies with such extensive spectral coverage. 
GM Aur is a K5 source in the Taurus molecular cloud, initially suspected to be a transitional disk 
(Marsh fc Mahonevl Il992l : iKoerner et al.l 1 19931 : iGhiang fc Goldreichl Il999l ) and later established as 
having an inner clearing of optically thick dust, though with a small amount of submicron- sized 


optically thin dust remaini ng, using Spitzer Infrared Spectrograph (IRS) observations (|Calvet et al 


20051 : lEspaillat et al.ll201()l i. Sub-milli meter observations confirmed that GM Au r has an inner hole 
in the dust, with a radius of ~ 20 AU ( Hughes et al. 2009 : Andrews et al. 2011 1. Rice et al. ( 2003 1 
showed that transitional disk IRS spectra may be produced by clearing of disk material by planets 
over a range of disk radii. While photoevaporation models can also produce inner disk holes, the size 
of the cavity in GM Aur and the p resenc e of remaining gas and small dust are better understood in a 
planet formation scenario (Owen et al. 2010|). In addition to existing long wavelength observations 
of GM Aur, far and near ultraviolet (FUV and NUV, respectively) spectra have provided glimpses 
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of the gas in the disk of GM Aur (lArdila et al. 2013|) and the emission produced in an accretion 
shock at the stellar surface (jlngjeby et al. 2013)- 


In g2l we discuss the UV, optical and near- IR observations used in our analysis. In g3] 
we outline our modeling procedures for each epoch of GM Aur data, including accretion shock 
modeling, comparison of IR fluxes to models of the emission from dust in the disk, and a discussion 
of the molecular contribution to the FUV. Finally, in g^l we propose a scenario to explain the 
variability of these spectra. 


2. Observations 
2.1. STIS 

GM Aur was observed during 3 epochs, referred to as GM_1, GM_2 and GM_3 in the following 
analysis, between September 2011 and January 2012 in GO program 11608 (PI: N. Calvet) with 
the Space Telescope Imaging Spectrograph (STIS). The observations were separated by 1 week 
(GM_1 to GM_2) and 3 months (GM_1 to GM_3) in order to study variability on short and long 
timescales. Spectra were obtained from the FUV to optical, approximately 1100 A to 1 /im, utilizing 
a combination of STIS gratings. The UV gratings, G140L and G230L used the MAMA detector 
and provided spectra between 1150 and 3180 A with resolution, R , ranging between 500 and 1,440. 
G430L and G750L used the CCD detector and covered 2900 to 10,270 A with R ~ 530 — 1, 040. We 
used data products produced through the STScI calstis reduction pipeline for our analysis. A log 
of all STIS observations is provided in Table Q] and the spectra are plotted in Figure [0 combined 
with IR data obtained from SpeX, discussed in H2.2I 


2.2. SpeX 


GM Aur was observed using the SpeX spectrograph (Rayner et al. 2003) on NASA’s Infrared 
Telescope Facility (IRTF). The SpeX observations were obtained using the cross-dispersed (here¬ 
after XD) echelle gratings in both short-wavelength mode (SXD) covering 0.8-2.4 /jm with a reso¬ 
lution of ~2000 and long-wavelength mode (LXD) covering 2.3-5.4 [in i and resolution ~2500. The 
LXD data were obtained on 5 nights (2011 September 11, 18, and 29 (UT) and 2012 January 2 and 
6 (UT). On three of these nights (November 11 &; 18 and January 6) SXD data were also obtained. 
For all XD observations, a 0.8" wide slit was used. The spectra were corrected for telluri c ext inc¬ 


tion, and flux calibrated a gainst the A0V star HD 25152 using the Spextool software (jVacca et al 


20031 : [Cushing et ah 20Q4) running under IDL. On nights when both SXD and LXD spectra were 
obtained, the LXD spectrum was merged to the SXD, by re-scaling in the LXD spectra in the 
wavelength region where they overlapped. 


Due to the light loss introduced by the 0.8" slit used to obtain the XD spectra, changes in 
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Fig. 1.— UV, optical and IR spectra of GM Aur. Three epochs of observations are shown as solid 
black (GM_1), dashed blue (GM_2) and dot-dashed red (GM_3) lines. Gaps in the near- IR SpeX 
data are regions where data is not shown due to poor telluric subtraction. 
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telescope tracking and seeing between the observations of GM Aur and a calibration star may result 
in merged XD spectrum with a net zero-point shift compared to their true absolute flux values. To 
correct for any zero-point shift that results from the 0.8" echelle observations, we also observed GM 
Aur on 4 of the nights (excluding only 2011 September 29, due to poor sky conditions) with the 
low dispersion prism in SpeX using a 3.0" wide slit. On nights where the seeing is 1.0" or better, 
this technique yields absolute fluxes that agree with aperture photometry to within 5% or better. 
The resulting scaling factors derived for the XD spectra were 10% or less on all 4 nights, and the 
scalings were applied to each spectrum prior to analysis. The three epochs closest in time to the 
STIS spectra were included in our analysis and details of the SpeX observations are listed in Table 

m 


3. Analysis and Results 
3.1. Accretion rates from STIS spectra 

It is generally accepted that the NUV and blue optical excess emission in classical T Tauri 
stars (CT TS) is produced in the a c cretion shock a t the stellar sur face, a result of magnetospheric 


accretion (jUchida &; Shibatal 1 19851 : iKoenigll Il99ll : IShu et al.1 119941 ) . In this scenario the circum 


stellar disk is truncated by strong stellar magnetic fields at a few stellar radii, inside of which 


material is channeled along the field lines to the star. Calvet & Gullbring (J1998) described the 
emission from an accreting column of material which forms a shock at the stellar surface. The 
shock emits soft X-rays towards the star, heating the photosphere below, and along the column of 
accreting material where it is reprocessed producing a spectrum peaked in the NUV. Models which 


(Herczeg & Hillenbrand 

2008; 

Rigliaco et ah 

2011- 

Manara et ah 


2014 ). As a second o rder approximation to the accretion shock model of lCalvet &; Gullbrind (119981 ) . 
Ingleby et ah (2013j) allowed columns with a variety of energy fluxes, J-, and surface coverage or 
filling factors, /, to co-exist. The presence of multiple accretion columns was motivated by models 
of the magnetospheric geometry which include higher order moments of the ste l lar magnetic field - 
producing a complex ac cretion pattern at the stellar surface (jDonati et al.l 120081 : iLong et al.l 120081 : 
Gregory fc Donati 2011 ). 


We fit the STIS NUV and optical spectra with the multi-column accretion shock model 
described in llnglebv et al.l (120131 ). STIS spectra were first de -reddened using Ay = 0.8 from 


Espaillat et al.l ( 201ll ) and the extinction law towards HD 29647 ( Whittet et al.ll2004J ) . We account 


for the emission from a T Tauri chromosphere and photosphere by using a non-accreting T Tauri 
star oi_weak T Tauri star (WTTS) template, specifically a STIS spectrum of the WTTS RECX 1 
from llnglebv et ah ( 2011 ). RECX 1 is in the nearby r/ Chamaeleon star forming region and has the 
same spectral type as GM Aur (K5). 


To scale the spectrum of the WTTS RECX 1 to that of GM Aur we use the veiling value at the 
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V band. Veiling refers to the degree by whi ch p ho tosp heric absorption lines appear shallower due 
to the addition of a continuum excess (Hartigan et ajj 1991), which we assume is accretion shock 
emission, or Fv^wtts = Fv,CTTs/( 1+^v), where Fv^wtts and Fv,ctts are continuum fluxes for the 
WTTS and CTTS, respectively, and ry is the amount of veiling, defined as ry = Fyy e u/Fy^wTTS, 
where Fy y e n is the flux of the veilin g co ntinuum. Since we lack high resolution observations, we 
use the value ry = 0.2, from lEdwards et al.l (|2006l ). Veiling is variable an d depends on th e amou nt 
of excess emission produced in the shock, therefore assuming the value of Edwards et ah (2006) is 
a source of error. However, we expect this error to be small given that the low accretion rate of 
GM Aur produces a small amount of veiling. The optical spectrum of GM Aur is dominated by 
the stellar emission, with the shock emission becoming important in the UV. 


As in llnglebv et al.l (2013), each accretion column is characterized by an energy flux, F, and a 
filling factor, /. The energy flux, T = 1/2 pv^, is a measure of the density of material in the accretion 
column, p, assuming that the magnetospheric truncation radius is not changing and therefore the 
infall velocity, v s is constant. Higher T columns deposit more energ y, pro ducing hotter shock 
emission and shifting the shock spectrum towards shorter wavelengths (Ingleby et ah 2013). The 
filling factor / gives the fraction of the stellar surface covered by the column. The stellar mass 
and radius are additional input parameters to the models. We estimate stellar properties using 
a synthetic J band magnitude found by integrating the 2MASS J band transmission curve over 
the de-reddened SpeX o bservation s. Th e sy nthetic J magnitude combined with the conversion to 
bolometric magnitude of lKenvon fc Hartmann! (119951 ) provides the stellar lu minosity, L t = 0,9 L^. 
We assume a distance to GM Aur of 140 pc, and a temperature of 4350 K (Kenyon & Hartmann 
199511 to g e t Il t = 1.7 Rq. Finally we estimate a mass M* = 1.1 M 0 using the isochrones of 


Siess et ah 


( 2000 ) 


We calculated models of the emission from accretion columns with log F= 10, 10.5, 11, 11.5 
and 12 erg s -1 cm -2 , corresponding to heated photosphere temperatures between 4600 and 10,700 
K. As F increases, the shock spectrum gets bluer, theref ore, a c ombination of columns provides 
the best fit to red and blue wavelengths. Following llnglebv et al.l (2013), we found the best model 
of the STIS spectra by calculating the reduced chi squared fit ( Xred ) °f combinations of accretion 
columns, each scaled by /, left as a free parameter. The accretion shock models do not attempt 
to reproduce line emission, therefore, the y 2 erf fit omitted spectral regions with strong emission 
lines. The assessment of the fit also excluded the region between 3200 and 3400 A, where there is 
significant noise due to edges of the NUV and optical gratings. To produce the final fit, we added 
the emission from each contributing accretion column to the emission from the WTTS template, 
RECX 1. Figure [2] shows the best fits to each epoch of STIS observations. 


The filling factor of the contributing columns at each epoch is given in Table [2j We list a 
characteristic value for the energy flux, F w , which is an average of the energy flux of all columns, 
weighted by the filling factor of each column. The listed total filling factor, ftotal, is the sum of 
the / values for the contributing columns, providing the total percentage of the surface covered by 
accretion columns. 
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Fig. 2.— Accretion shock model fits to 3 epochs of GM Aur STIS spectra. In each panel the black 
line is the de-reddened STIS spectrum, while the blue line is a STIS WTTS template, here RECX 
1. The broken lines represent accretion columns characterized by different values of log J-= 10, 
10.5 (short dashed), 11 (dot-dashed), 11.5 (dash-triple dotted), and 12 (long dashed) erg s _1 cm -2 . 
No columns of log J ~=10 erg s -1 cm -2 were needed for the fits. The red line is the best fit to the 
data of the WTTS template plus the summed accretion column emission. The gray regions in the 
first panel represent spectral regions excluded when calculating the model fit to the data. 

































The mass accretion rate was found by by adding the contributions of the columns (Ingleby et al 


2013) 


M = 


8ttRz 


X^i/i 


8vr Rl ^ t 

9 wj total 


(i) 


We assume that v s is 0.89 of the free-fall velocity, which for the stellar parameters is 494 kms -1 , 
corresponding to infall from ~ 5 stellar radii. Values for M in the three epochs are given in Table El 
The error in our estimates of the accretion rate is dominated by uncertainti es in A y. Existing esti¬ 


mate s of extinction for GM Aur range from 0.1 (jKenvon &; Hartmannlll995l) to 0.8 ( Espaillat et al 


2011) with a value of 0.3 derived fro m long waveleng t h opt ical spectra (IHerczee &; Hillenbrand 


2014). We adopted the value listed in lE s pailla t et aL (2011) in order to directly compare those 


results to our analysis. The range of measured Ay values introduces a factor of 2 error in our M 
estimate. 


3.2. FUV to IR variability 

Here, we analyze the wavelength dependent variability in the STIS and SpeX data and dis¬ 
cuss some of the common contributors to the emission in each spectral region. GM_1 and GM_2, 
separated by one week, have very similar spectra and in TS.ll we showed that these two epochs had 
minimal differences in the accretion properties. However, GM_3, which was observed 3 months after 
GM_1, is clearly dimmer, though the magnitude of the dimming across the spectrum is not uniform 
(see Figure Q]). Namely, the biggest decrease in fluxes occurs in the NUV, with accompanying yet 
smaller dips in FUV, optical and IR fluxes. In the following analysis we focus on the differences 
between GM_1 and GM_3. 


3.2.1. NUV and optical excess 

As shown in Table El while the total filling factors for all three epochs are approximately equal, 
the energy flux of the columns is higher in GM_1 and GM_2 than in GM_3. As a result, the columns 
in GM_3 emit less in the NUV, significantly lowering the total flux in this wavelength range. This 
can best be seen in Figure [3l which compares the normalized total flux due to the accretion columns 
in the three epochs. 

Figure |4] shows the difference in flux between GM_1 and GM_3, with excesses observed in GM_1 
at all wavelengths. We fit the NUV and optical portion of this difference spectrum with the emission 
from an accretion column. We tried fits with columns characterized by log T = 10 — 12 and found 
that the best fit to the excess comes from an accretion column with T = 3 x 10 11 erg s^ 1 cm~ 2 
and / of 0.006. The accretion rate of a column with these parameters is 5.1xl0~ 9 -ATq yr . This 
M is roughly consistent with the difference in M between GM_1 and GM_3 listed in Table El 

To summarize, we find that the difference between the accretion properties of GM_1 and GM_3 
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Fig. 3.— Total shock models used to fit the STIS spectra. The contribution from each accretion 
column in the fit is summed for the final accretion shock models. Fluxes are scaled for direct com¬ 
parison of the shape of the shock emission. The emission from the accretion shock is characterized 
by higher density columns for GM_1 (solid) and GM_2 (dashed), or higher T. GM_3 (dash-dotted) 
has the lowest density of accreting material and therefore the coolest shock spectrum. 
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Fig. 4.— Excess in GM_1 over GM_3. The spectrum of GM_3 was subtracted from GM_1 to isolate 
the excess emission observed in the first two epochs (solid black line). The excess is fit in the 
NUV and optical with accretion shock emission characterized by T = 3 x 10 11 erg s _1 cm -2 and 
/ = 0.006 (dashed purple line). The optically thin dust emission is also shown to contribute to 
the IR emission (dash dotted, magenta line). The solid red line is the total of the accretion shock 
emission and the optically thin dust emission. Here we are not attempting to explain the FUV 


excess. 
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is in the energy flux, or density, of the accretion columns rather than the total surface coverage. 
The decrease of M in GM_3 is mainly due to the decreased density in the high density accretion 
columns. 


3.2.2. IR excess 


GM Aur is known to be variable in the infrared. Using Spitzer IRS spectra, lEspaillat et al 


(2011) found a change in the strength of the 10 pm silicate emission and in infrared continuum 
emission from 5 to 10 pm. They interpreted this as due to changes in the amount of optically 
thin dust within the disk hole. Espaillat et ah (2011) found that the inner disk hole contained 
2xl0~ 12 Mq of optically thin dust and they could reproduce the observed variability by decreasing 
the amount of dust in the optically thin region by a factor of two. 

In Figure |4] we detect an excess in the IR in GM_1 as compared to GM_3. This excess is 


(Fischer et ah 

2011: 

McClure et ah 

2013) 


to fit the SpeX data for GM_1 and GM_3, following the procedure of Espaillat et al. (2010, 2011) 
to reproduce the excess. The dust is heated by stellar radiation and we change the mass in dust 
to fit th e spec tra. The silicate dust composition used in the optically thin region is the same as 


m 


Espailla t et ah (2Q11), which is based on the best fit to the 10 pm silicate feature. We use 


ISM-sized (0.005-0.25 pm) dust composed of 32% organics, 12% troilite, and 56% silicates. The 
silicates are 90% amorphous and 10% crystalline forsterite. We find that to fit the GM_1 spectra 
we need 1.9xl0 -12 M 0 of dust in the optically thin region, while we need 1.4xl0 -12 M 0 for GM_3. 
Figure [4] compares the GM_1-GM_3 excess spectrum to the emission from the difference between 
the two optically thin dust models used to fit the SpeX data. 


3.2.3. FUV excess 


Figured] shows that the observed FUV excess in GM_1 compared to GM_3 cannot be explained 
by emission from more energetic accretion columns alone. However, it is well known that in ad¬ 
dition to shock emission, molec ular gas emission contributes to the FUV spectrum, primarily Ho 
with a contribution from CO ( Herczee et ah 2002 . 2004 : Bergin et ah 2004 : Inglehv et ah 2009 : 
France et al.ll2011al lb[h This molecu lar gas was expected to be within a few AU of th e central star, 


as Ho was spatially unresolved by iHerczeg et al 


((2003) for TW Hya, and indeed llnglebv et al 


(2011) showed that spectrally resolved H 2 lines originated between the magnetospheric truncation 
radius, out to ~3 AU for the case of RECX 11 in p Chameleon. 

H 2 in the inner disk is excited via one of two mechanisms. The first is excitation by Lya 
photons populating the upper levels of H 2 which then creates a fluore scen t spectrum as it de- 
excites. While Lya fluorescent lines were identified by IHerczeg et al. (2002) in high resolution 
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STIS spectra, the features are detectable at low resolution and are found in the GM Aur FUV 
spectra for all epochs. The second mechanism is collisional H 2 excitation by free hot electrons. 
Metals present in the inner disk are ionized by X-rays from the star and these electrons ionize the 
abundant H and He producing a large supply of hot electrons. Electrons collisionally excite H 2 and 
therefore different paths for de-excitation are allowed with one path producing continuum emission 
as the H 2 is dissociated. Therefore the two excitation methods produce unique and distinguishable 
H 2 emission spectra (lAbgrall et al.1 1997). 
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Fig. 5.— FUV excess in GM_1 - GM_3 after shock subtraction. The black line shows FUV 
differe nce spec trum minus the emission from the accretion column. A continuum feature discussed 


m 


Bergin et ah (.2004, B04) is identified near 1600 A and produced by electron impact excitation 


of H 2 . The wavelengths of several lines produced by Lya fluorescence are marked by red dashes. 

3 


In Figure [5] we show the excess FUV emission after the subtraction of the accretion shock 
continuum emission (black minus red in Figure [4]). The strong FUV atomic lines remain after the 
removal of the shock, since the accretion shock models do not attempt to produce line emission; 
however lines such as the Si IV doublet near 1400 A, C IV at 1548 and 1551 A and He II A1640 
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A are known to co rrelate with the mass accretion rate ( Johns-Krull et al.ll2000l : ICalvet et al .1120041 : 
Ardila et al.ll2013 ). In addition to the strong hot lines, we see a remaining continuum emission which 
cannot be explained by the accretion shock. The broad continu um feature at 1600 A was identified 


as col l isiona lly excited H 2 emission from the disk of CTTS in iBergin et al.1 (120041 ). Ilnglebv et al 


(2009, 2012) showed that this feature is always observed in CTTS but absent in WTTS, when 
the disk has been dispersed. We observe this continuum feature in the FUV difference spectrum, 
revealing an excess of collisionally excited H 2 emission coming from the disk in GM_1. Assuming 
that the difference reflects an increase in H 2 density, we estimate the necessary density difference 
assuming that the flux in the 1600 A H 2 feature (Figoo) is proportional to the H 2 surface density 
(X# 2 ), or T\eoo °c (Ilngleb v et a l. 20091 ). Using this relation, we find that the surface density 
of H 2 decreased from GM_1 to GM_3 by a factor of 1.4. 

The difference spectrum does not support an excess in the Lya excited H 2 emission. Flu¬ 
orescent lines easily detected in the original spectra do not appear in Figure [5l indicating that 
there is little change in this type of H 2 emission when the accretion rate drops. In addition to the 
observations discussed here, GM Aur was observed with STIS in April, 2003 in GO Program 9374 
(PI: E. Bergin) and during August, 2010 in GO Program 11616 (PI: G. Herczeg). The 2010 spectra 
were analyzed and accretion rates determined in Ilnglebv et al.l (12013! ) while we calculate M for the 
2003 data using the method described in 43.11 Accretion rates in 2003 and 2010 were 9.1 x 10~ 9 
and 9.6 x 10~ 9 A/gyr -1 , respectively. We compared the fluxes in the strong Lya fluoresced H 2 
lines of all 5 STIS epochs and found that there is no evidence for a correlation with M, based on a 
Pearson r coefficient between 0 and 0.5 for most lines. While Lya is heavily absorbed by interstellar 
hydrogen and cannot be observed directly, these results suggest Lya is not changing between GM_1 
and GM_3, contrary to predictions that Lya is formed in the accretion shock. If instead Lya is 
formed in the broad accretion flows, then perhaps the magnetospheric truncation radius is larger 
when M drops, increasing the Lya emitting region. While interesting, much more substantial work 
is needed to interpret Lya in CTTS, for example with efforts to reconst ruct the Lya line profile by 
observing H 2 fluorescence (IHerczeg et al.ll2002l : ISchindhelm et al.ll2012l L 


4. Discussion 

In the previous sections, we presented observations which displayed variability in the FUV-NIR 
emission of GM Aur. From these observations, we extracted four measurements of the properties of 
GM Aur. First, we measured a lower accretion rate in GM_3 than in GM_1 and GM_2. For GM_3, 
we also measured a lower dust mass and a lower H 2 density in the inner disk than the earlier GM_1 
and GM_2 data. Lastly, the change in the FUV-NIR emission of GM Aur was observed to occur 
over a timescale of 3 months. In the the following, we consider which physical properties of the 
system may have been different during GM_3 to explain our observations. In particular, we focus 
on how inhomogeneities in the inner disk may explain the four measured properties of GM Aur. 

We can explain the decrease in the NUV emission of GM Aur with our accretion shock models, 
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which measure a lower accretion rate in GM_3. It is important to define what a single measured 
accretion rate corresponds to. Accretion rates are derived by measuring the shock emission in 
excess above the stellar photosphere and using the accretion shock model (or a slab model) to 
account for the emission outside the observed wavelengths (i.e., <1000 and > 1 fin i in this work). 
Determining the accretion luminosity from t he excess emission, a n d using the stellar mass and 


radius, the accretio n rate then follows (e.g., iGullbring et al.1 Il998l : iHerczeg fc Hillenbrand 


2008 


Inglebv et al.l 2013|). And so, the reported accretion rate is a measure of the emission created 
from mass impacting a region of the stellar surface at a particular point in time, which in turn 
corresponds to a particular orientation of the star (and hence magnetosphere) and disk along our 
line of sight. Analysis of the spectral energy distribution of the excess emission in 113.11 indicates 
that in GM_3 we were observing a magnetospheric region characterized by low density accretion 
columns. This does not preclude that high density columns were present at other azimuthal regions 
that were not visible at that time. 

Considering 5 epochs of existing STIS spectra (combining the current observations with those 
obtained in 2003 and 2010), 4 out of 5 times GM Aur was observed during a high state of accretion 
and in GM_3 we caught it in a low state. Compared to results from accretion shock fitting of 
our multi-epoch data, the accretion rates in 2003 and 2010 were similar to the higher accretion 
rates of GM_1 and GM_2 (see 113.2.31) . Specifically, in the multi-column fitting of the spectra, all 
but the GM_3 ep och are c hara cterize d by high density accretion columns with small filling factors. 
As discussed in I nglebv et al.l (120131 ). lower density columns could still be contributing, but their 
emission could not be extracted from the brighter stellar spectrum. Whereas during GM_3 the 
accretion columns had lower densities with larger stellar surface coverage. Although schematic, 
this evidence points to an inhomogeneous and variable magnetosphere. This is not surprising given 
the highly complex natur e of the stellar magnetic fields d etermined by Zeeman Doppler techniques 
(jDonati et al.ll2008l . l201ll . 120131 : iGregorv fc Donatill201ll ). However, in addition to the complexity 
of the magnetic fields, the structure of the inner disk from which the mass is coming from must 
play a role in the driving the full spectrum of variability. 

Simultaneous with the NUV decrease in GM_3, there are decreases in the NIR spectrum and 
FUV H 2 emission feature at 1600 A. We attribute these to lower dust masses and lower H 2 densities 
in the inner disk, respectively. In 113.2.21 we discussed the optically thin dust models used to fit the 
NIR excess and found that the dust mass was a factor of 1.4 lower in GM_3 than in GM_1. Since the 
dust is optically thin, this means the dust surface density decreased by a factor of 1.4. Interestingly, 
1.4 is the same factor by which we estimate the H 2 density decreased between GM_1 and GM_3 in 
113.2.31 This suggests that the H 2 1600 A feature is tracing the same location in the disk as the dust, 
and that the dust and gas are well-mixed in the inner disk. The decrease in the H 2 emission feature 
at 1600 A is not accompanied by a change in H 2 fluorescent lines, which should also trace the gas; 
however, this may be due to different origins of the H 2 lines. As discussed earlier f 113.2.311 . the H 2 
emission feature at 1600 A is a c onsequence of collisional excitation fo llowing X-ray ionization of 
heavy elements in the inner disk (jAbgrall et al.lll997l : iBergin et al .1120041 ) while fluorescent lines are 
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due to excitation of electrons to high electronic levels of H 2 by UV Lya emission (IHerczeg et al 
2002 1. X-rays penetrate deeper into the disk than UV photons because U V rad iation is absorbed 
by small dust grains in the uppermost layers of the disk (Nomura et ah 20071 1. We expect that 
the H 2 emission feature at 1600 A is coming from deeper in the disk and is tracing the same disk 
regions as seen in the NIR dust emission. 

Our analysis suggests a link between the mass in gas and dust in the disk and the mass being 
loaded onto the star at the magnetospheric truncation radius. However, the kinematic timescales 
where the gas and dust reside in the disk are problematic given the 3 month time difference between 
GM_1 and GM_3. At_the_dust des truction radius, the viscous timescale is ~ 2-10 years using 
Equation 37 of Hartmann et al. (1998) with Ri = 10 R*, T=5000-1000, and a=0.01. It is unlikely 
that the mass in the entire inner disk is decreasing within 3 months and that this is causing the 
lower accretion rate in GM_3 given that the dust and gas in the inner disk are fed by the outer 
disk on the viscous timescale. A period of 3 months requires a change in the disk to occur near 
the magnetospheric radius, ~5 R *. At this radius, the viscous timescale is no longer relevant as 
material is transported inward along magnetic field lines as opposed to viscously through the disk. 

In order for a similar decline of mass in the dust and H 2 , the two must be spatially related, 
but dust must reside outside the dust destruction radius (typically ~ 10 R*) where the viscous 
timescale is much longer than 3 months. However, if we allow large dust grains, 1 mm, and a 
sublimation temperature of 1600 K (as opposed to the co mmonly used 1400 K), dust grains can 
survive down to 0.07 AU, fo llowing iD’Alessio et al.1 (120051 1. which is approximately the corotation 
radius. Stauffer et al. (2015) found that, for a sample of young stars, IR variability indicates that 
inner disk walls may be present near the corotat ion radi us. I n our simple model, we assume the 
magnetospheric truncation radius is 5 R* (I Cal vet & Gullbring 1998! ). or 0.04 AU. In reality, the 
truncation radius depends on the strength of the magnetic field and accretion rate. In addition, 
magnetosphere s are co mplex, therefore the truncation radius varies and may extend to corotation 
at times flD onati e t al. 2011). If some of the NIR emission is produced in the large grains near the 
midplane (allowing for small grains remaining in the upper layers), then the gas and dust may be 


co-spa tial near the corotation radius and also linked to an extended magnetosphere. Blinova et al 


(2Q15) showed that if the magnetospheric truncation radius ( r m ) is near the corotation radius (r cor ), 
r cor ~ 1.4 r m , then accretion is chaotic and variable, whereas for r cor > 1.4 r m , accretion is stable. 
Unstable accretion occurs through short lived “tongues”, contributing to inhomogeneous accretion. 

While the NUV emission is related to variable accretion, FUV and NIR variability is rooted 
in the disk. We cannot conclusively say what is leading to inhomogeneities in the inner disk 
with our low cadence data set. One theory is that disk material builds up near the corotation 
radius before emptying quickly onto the star (D’Angelo & Spruit 2012 ). This would produce a 
period of lower density in the inner disk after the material is dumped onto the star. While this 
phenomenon is typ ically used to describe objects undergoing large accretion bursts (like EXors), 
D’Angelo fe Spruitl (2012) note that any object with a strong magnetic field and low accretion 
rate should undergo this affect at some level. As discussed, the magnetosphere may truncate the 
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disk closer to corotation than typically assumed, and gas and large dust grains may coexist near 
corotation, linking the decline in disk mass observed with the lower accretion rate. Alt ernat ively, 


i t is p ossible that inner disk gas may have an emission component in the NIR. Tannirkulam et al 


(2003) required the presence of NIR emission inside the dust destruction radius of Herbig Ae stars 
and suggested optically thick gas as the source. If the NIR does come from inner disk gas, then 
the inhomogeneities and magnetospheric truncation radius may reside inside corotation. 

Another common cause of variability in CTTS stems from variable extinction by dust. GM Aur 
becomes redder as it dims during the GM_3 epoch and objects with that behavior m ay be unde rgoing 
such variable dus t extinction, possibly from warps in the circumstellar disk (jGiinther et al.l 120141 : 
Codv et alJl2014l ). Variable obsc uration requires that the disk be observed close to edge on, as in 
the well st udied case of AA T au (jBouvier et al.l 120131 1. GM Aur has a well determined inclination 
of ~ 57° (Dutrev et al. 1998}), therefore dust obscuration requires disk structures at large scale 
heights, unlikely to exist for any sustained period of time. A dimming event in the CTTS RW 
Aur A (inclined between 45° and 60°) may be due to d ust ca rried to large scale heights above the 
disk in an enhanced wind, producing the obscuration (Petrov et al. 2015); however, RW Aur A is 
known to ha ve st rong outflows as e videnced by existing jets, whereas GM Aur shows no sign of jets. 
In addition, Elle r br o ek et ah (20141) showed that dust lifted to large scale heights above the disk, 
while producing dimming in the optical by obscuring the star, would increase the NIR emission as 
the surface area of dust is increased. 


In conclusion, we propose that the observations presented here are tracing inhomogeneities in 
the inner disk and magnetosphere. Parsing out the detailed structure of these reg ions and how 


thev interact is an extremelv difficult task given the high complexitv (e.g., ' 

domanova et al. 
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2013). In fact, there 


are some observed properties of GM Aur that we cannot yet explain. We note that while the dust 
mass and H 2 density decrease by a factor of 1.4 in GM_3, the accretion rate decreases by a factor 
of 2.8 (although given the errors, these numbers may in fact be more similar). The surface density 
may not scale exactly with M as factors such as temperature and viscosity also play a role. It is 
puzzling what causes this difference, given that it is not clear exactly how accretion occurs onto 
the star. In any event, our data connect the disk, near the footprints of the magnetic field lines, to 
the material in the accretion shock at the stellar surface. Higher cadence observations of GM Aur 
would be ideal to test our proposal and to make further progress on the issue of variability of TTS. 


5. Summary 

We analyzed variability in the FUV to NIR spectrum of GM Aur over 3 epochs of observation, 
separated by one week and 3 months. We found that GM_1 and GM_2 were nearly identical in 
fluxes. GM_3 was observed during an epoch of low fluxes across the entire spectrum. We drew the 
following conclusions: 
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• The accretion rate of GM_3 was lower than the first two epochs by a factor of 2.8. While the 
surface coverage of the accretion columns remained approximately constant, the energy flux 
in the columns decreased in GM_3. Energy flux is related to the density, so GM_3 had low 
densities in the accretion columns compared to the earlier epochs. 

• We fit SpeX observations, obtained contemporaneously with the STIS spectra, with models 
of optically thin dust in the inner cavity of GM Aur. During the epoch of low accretion, the 
mass of dust in the disk decreased by a factor of 1.4. 

• The FUV spectrum indicates that the emission from H 2 molecules decreased during the epoch 
of low accretion and decreased dust mass. In particular, we found that the spectrum of H 2 
collisionally excited by free electrons in the disk described the shape of the excess in GM_1 
over GM_3. Moreover, we find that the mass of H 2 changed by the same factor as the dust. 

• GM Aur is seen here to vary within a 3 month period, placing the origin of variability near 
the magnetospheric truncation radius. If the disk is truncated near corotation and dust can 
survive near corotation in large grains, then inhomogeneities at that location in the disk can 
account for accretion variability. The build up and subsequent dumping of material onto the 
star is one possible explanation for the source of disk inhomogeneities. 

We conclude that these FUV-NIR data are providing a glimpse at the complexity of the 
interaction between the stellar magnetic field and the inner disk. Future observational work is 
needed to test this conclusion as well as theoretical work on the observable properties of star-disk 
interactions. 
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Table 1. Log of Observations 


Epoch 

Telescope/ Instrument 

Date of Obs 

GM_1 

HST/STIS 

09-11-2011 

GM_1 

IRTF/SpeX 

09-11-2011 

GM_2 

HST/STIS 

09-17-2011 

GM_2 

IRTF/SpeX 

09-18-2011 

GM_3 

HST/STIS 

01-05-2012 

GM_3 

IRTF/SpeX 

01-06-2012 



Table 2. 


Results from Multi-Component Model Fits to UV and Optical Spectra 


Object 

/(io 10 ) 

/( 3 X IO 10 ) 

/(io 11 ) 

/( 3 x 10 11 ) 

/(io 12 ) 

ftotal 

Weighted T 

M (M 0 yr x ) 

GM_1 

0 

0 

0.03 

0 

0.001 

0.031 

1.3 x 10 11 

1.1 x IO" 8 

GM_2 

0 

0 

0.02 

0 

0.001 

0.021 

1.4 x 10 11 

8.5 x IO" 9 

GM_3 

0 

0.04 

0 

0.0006 

0 

0.041 

3.4 x 10 10 

3.9 x IO" 9 


Note. — The values in parentheses represents the energy flux (J 7 ) of each column in erg s 1 cm 








